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ABSTRACT 



We present an analysis of the frequency shift and the even 
terms of the frequency splitting coefficients carried out 
using six years of LOWL data, starting in 1994. The 
temporal variations, and their dependences with the fre- 
quency and degree are addressed. The results are consis- 
tent with previous analysis. 



1. INTRODUCTION 



It is very well established that the variations of some of 
the oscillation mode parameters is a signature of the solar 
activity cycle. For instance, the acoustic mode frequency 
is shifted positively of 0.4//Hz from minimum to max- 
imum of the solar cycle (Libbrecht and Woodard 1990, 
Regulo et al. 1994, Jimenez-Reyes et al. 1998). It has 
also been noted that, this variation depends on the fre- 
quency and on the degree of the mode as well. The even 
terms of the splitting coefficients present similar changes. 
Here, we present an analysis of the frequency shift using 
LOWL observations. 



2. DATA ANALYSIS AND RESULTS 



The raw data consist in full disk solar Doppler images, 
which have been collected by LOWL instrument. This 
experiment is based on a Magneto-Optical Filter and has 
demonstrated its high sensitivity to the solar oscillations. 
The data were analyzed using the LOWL pipeline which 
has recently been improved (Jimenez-Reyes 2000). The 
images are first calibrated and then a spherical harmonic 
decomposition is performed in order to create time series 
for degrees from £=Q to 99. Finally, a Fast Fourier Trans- 
form is applied to the time series. 

The spherical harmonics are not orthogonal over the ob- 
served area, limited to one hemisphere. Therefore, the 
modes cannot be totally isolated and some correlations 
between different (I, m) spectra exist. Historically, the 



X 2 with 2 degrees of freedom. However, this assump- 
tion is not correct when the spectra are not independent. 
Recently, important improvements have been achieved in 
the fitting techniques to be applied to observations made 
using instruments with spatial resolution (Schou 1992, 
Appourchaux et al. 1998). The real and the imaginary 
part of the Fourier Transform follow a multi-normal dis- 
tribution which is described by a covariance matrix. The 
likelihood function under these assumptions can be writ- 
ten as: 
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where yivi) corresponds to the Fourier Transform and 
E(a, Vi) is the covariance matrix calculated by: 
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where C is the leakage matrix which represents the cor- 
relation between modes with different m (m-leakage), 
whereas B is the noise covariance matrix which gives 
information about the noise correlation between spec- 
tra. v(a, v) is the variance and it is given by a simple 
Lorentzian profile defined by the parameters a that we 
are trying to infer using a maximum likelihood method. 
The shift between m-components is given by: 

Tlcoef 

Vnlm = Vnl + ^ a i( n > ( m ) ( 3 ) 

where Pf{m) are orthogonal polynomials normalized 
such that Pf(£) = £ (Schou et al. 1994, App. A). The odd 
a-coefficients are induced by the internal rotation while 
the even terms are mainly related to effects of second or- 
der due to rotation, the presence of magnetic fields or any 
departure from spherical symmetry. 

The distance between modes with the same order n and 
with A£=±l get closer and closer at high degrees lead- 
ing to a leakage of the modes (^-leakage). This prob- 
lem arises at lower degree in the case of observations 
from just one station, due to the existence of the side- 
lobes as a result of the modulation of one day in the sig- 
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v/L « 31/iHz. In order to reduce the systematic er- 
rors an extra covariance matrix is added to take into ac- 
count the presence of spurious modes which can be close 
to the target mode. Nevertheless, we have seen that the 
systematic errors are specially important in the odd a- 
coefficients and not in the even ones. 

At high degree, the numerical evaluation of the likelihood 
function get slow, mainly because of the large dimensions 
of the covariance matrix. Therefore, we use just the main 
diagonal of the covariance matrix. We have checked that 
this is indeed a very good approximation for high degrees. 



3. FREQUENCY SHIFT 
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Figure 1. Time variation of the integrated frequency shift. 
The BBSO data have been averaged in the same way and 
plotted in the figure. In addition, 15 years of Mark-I fre- 
quency shifts are shown as well. 

Figure [l] shows the temporal variation of the integrated 
frequency shift for frequencies between 2 and 3.5 mHz. 
For comparison, the figure includes the same calculations 
using the BBSO data (squares). Mark-I instrument has 
been measuring the solar oscillation for more than 20 
years. This instrument observes the Sun like a star, so 
the information is limited to low degrees (£ < 3). Re- 
cently, the database has been re-analyzed (Jimenez-Reyes 
2000) over a much wider time interval than before. Here, 
we have calculated the yearly frequency shift every six 
months for the last 15 years . All the results are very well 
correlated with the changes in the solar activity cycle de- 
noted on Fig. |l] by the best linear fit between the radio 
flux and the results from Mark-I. The amplitude of these 
changes are, for low degrees, close to 0.4/iHz, peak-to- 
peak, and the correlation with the radio flux variations is 
close to 1 . 

The temporal variations of the central frequency are ex- 
pected to change with both frequency and degree. Ear- 
lier works (Libbrecht and Woodard, 1990) have shown a 
strong variation with the frequency, the frequency shift 
being null at ^2 mHz and increasing at higher frequency. 
The ^-dependence has been found to be more weak. 

Figure || shows the variation of the frequency shift av- 
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Figure 2. Frequency dependence of the frequency shift for 
three years. The error bars have been shown just for one 
year for clarity. The solid line represents the best linear 
fit of the inverse mode mass to the 1999 observations. 



the mean over the six years of data. The last minimum 
of the solar cycle happened in 1996, whereas 1999 rep- 
resents the year with higher level of solar activity in our 
database. As we can see, the frequency dependence is 
very clear, mainly between 1996 and 1999. The inverse 
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Figure 3. t-dependence of the frequency shifts for three 
different years. The solid line through the results corre- 
sponding to 1999, denotes the best linear fit to the inverse 
mode mass. 

mode mass has been calculated using the model of Morel 
et al. (1997). It has been averaged in the same way as 
the results shown in Figure |2|. The solid line over the 
points corresponding to 1999 represents the best linear 
fit to inverse mode mass. The observed frequency shift 
exhibits almost the same frequency dependence than the 
inverse mode mass. However, there are still some differ- 
ences which can be seen also in the analysis of the GONG 
data (Howe et al. 1999). This may or not be significant 
but it is interesting to notice that both results present the 
same fluctuations around the inverse mode mass curve. 

The ^-dependence has been also calculated and plotted in 
Figure |[ Again, we have selected the same years than 
in the last figure and the reference is the mean over the 
six years. The dependence with the degree is maximum 
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and £=99. The inverse mode mass, averaged in the same 
way than the observations, is shown in the figure, scaled 
to the best linear fit using the results corresponding to 
1999. The solid line follows remarkably well the fre- 
quency shift. 

The variations of the frequency shift have been studied, 
assuming a linear relationship with the solar activity in- 
dices. 

The solar indices used here are: 

• the sunspot number,R/; 

• the integrated radio flux at 10.7cm, Fio 

(both obtained from the Solar Geophysical Data); 

• the Kitt Peak magnetic index (KPMI) extracted 
from the Kitt Peak full disk magnetograms (Harvey 
1984); 

• the Mount Wilson index also called Magnetic Plage 
Strength Index, MPSI (Ulrich et al. 1991); 

• the equivalent width of Hel 10830A averaged over 
the whole solar disk using data from Kitt Peak ob- 
servatory; 

• the Total Solar Irradiance, TSI (Frohlich, C. & Lean, 
J., 1998). 

The results shown in Tab. 1 can be compared with those 
obtained recently by Howe et al. (1999) and Bhatnagar 
et al. (1999) from GONG observations. The magnitude 
of these parameters are very similar although some dif- 
ferences remain. The reason could be the different range 
in frequency taken to calculate the average. While we 
are taking frequencies between 2 and 3.5 mHz, Bhatna- 
gar et al. (1999) were taking all the result between 1500 
and 3500/iHz and Howe et al. (1999) were using a mean 
calculated at 3mHz. 

Table 1. Results of the weighted linear least-squares fits 
for the frequency shifts as a function of different solar 



indices. 


Index 


Intercept (/iHz) 


Slope 


Rj 


-0.0314±0.0070 


0.0031±0.0001 a 


Fio 


-0.2245±0.0131 


0.0032±0.0001 b 


KPMI 


-0.2203±0.0191 


0.0317±0.0018 c 


MPSI 


-0.0168±0.0063 


0.1604±0.0068 d 


He 


-0.5286±0.0416 


0.0114±0.0008 e 


TSI 


-363.4905±52.6823 


0.2662±0.0386 / 



The results for two selected years, (1996 in the bottom 
and 1999 in the top) are shown in Figure The even 
a-coefficients are close to zero in 1996, when the solar 
cycle reached the last minimum. However, a big change 
can be seen in 1999, following again the same frequency 
dependence than the inverse mode mass. Again, there is 
a clear fluctuation of a-i around the inverse mode mass 
curve. Nevertheless this may not be significant as it stays 
within the error bars. 
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Figure 4. Frequency dependence of the even a- 
coefficients for two different years, 1996 (bottom) and 
1999 (top). 

The time-variation of the solar cycle does not affect 
only the central frequency. The changes in the even a- 
coefficients demonstrate that, whatever is the perturba- 
tion leading to these variations, it depends on the latitude. 

We have integrated the even terms for each one of the 
time series. The resulting values are shown in Figure [| 
against the radio flux calculated over the same period. 
The straight line represents the best linear fit of both 
quantities. 
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4. THE EVEN A-COEFFICIENTS 



The even terms are expected to change with the fre- 
quency. Therefore, they were ^-averaged over short re- 



5. CONCLUSION 

We have analyzed six years of LOWL data from 1994 to 
the end of 1999 in order to parameterize the variation of 
the frequency shift and the even a-coefficients. 

In summary, the acoustic central frequency presents a 
variation which is very well correlated with the solar cy- 
cle. The integrated amplitude of this variation, peak-to- 
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Figure 5. Variation of the integrated even a-coefficients against the radio flux at 10.7. cm. The solid line shows the best 
linear fit. 



function of both the frequency and the degree, and these 
dependences can be well fitted by the inverse mode mass. 

In addition, the even a-coefficients present a significant 
change with the solar cycle, very well correlated with the 
solar indices. Moreover, as for the frequency shift, these 
changes increase with the frequency. The results agree 
very well with earlier studies (Libbrecht and Woodard, 
1990) carried out at the beginning of the last solar cycle. 

LOWL instrument has been recently updated and a new 
instrument has been developed and installed at the Ob- 
servatorio del Teide (Tenerife), leading to a new network 
called Experiment for Coordinated Helioseismic Obser- 
vations (ECHO). ECHO is intended to continue observ- 
ing for a complete solar cycle allowing us to better track 
the origin of the solar activity cycle. 
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